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We report on an exploratory study of Nf = 2 dynamical domain wall fermions and the DBW2 gauge action at 
weak coupling. Details of improved simulation algorithms and preliminary results for the hadron spectrum and 
renormalised light and strange quark masses will be presented. 



1. INTRODUCTION 

In the quenched approximation Domain Wall 
Fermions |1I2| (DWF) have been found to be 
an extremely successful approach to simulating 
QCD on the lattice. Crucial to this success is 
the fact that, when working at weak couplings 
(a -1 ps 2GeV) and using improved gauge actions 
such as the DBW2 action, the degree of explicit 
chiral symmetry breaking is very small for prac- 
tically useful sizes of the fifth dimension (O(10)). 

Early simulations of dynamical DWF, per- 
formed at relatively coarse couplings, suggested 
that L s m O(100) would be needed before the 
degree of chiral symmetry breaking was small 
enough to be acceptable. Here we will report on 
the progress of a preliminary study |3] of Nf = 2 
dynamical DWF in which we both adopt the 
DBW2 gauge action, and move to weaker cou- 
pling, in an attempt to find a region of parame- 
ter space where dynamical DWF simulations are 
practical. 

2. SIMULATION PARAMETERS 

All the results that will be presented were gen- 
erated using the DWF action with L s — 12 and 
M5 = 1.8, and the DBW2 gauge action with 
P = 0.80 on 16 3 x 32 lattices. Using the HMC 
algorithm we have generated three separate evo- 
lutions for bare masses of m/ = 0.02, m,f = 0.03 
and rrif — 0.04. Table ^ summarises the total 
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number of trajectories collected so far, together 
with the acceptance. Each HMC trajectory is of 
length 0.5 in HMC time and is split up into 50 
leapfrog integration steps for mass of 0.02 and 
0.03, and 40 integration steps for 0.04. 

Table 1 



Evolution details 



arrif 


trajectories 


acceptance 


0.02 


4716 


78 % 


0.03 


4785 


78 % 


0.04 


3445 


68 % 



3. ALGORITHMIC DETAILS 

The number of degrees of freedom of DWF 
grows with L s , but the number of physical de- 
grees of freedom does not. To cancel off this bulk 
divergence a set of Pauli-Villars fields, 

$t£>t( TO/ = l)D(m } = 1)$ , (1) 

is added to the DWF Lagrangian for dynamical 
simulations. Previous work has used two sets of 
pseudo-fermion fields to represent the DWF ac- 
tion: one for the fermion piece of the action and 
one for Pauli-Villars. The cancellation between 
these two terms is therefore only apparent after 
the average over the pseudo-fermion fields. Here 
we have used the fact that 

dct(£>t(i)j>(i)) 
det(Dt(m/)Z»(m/)) 1 ' 

is equal to 

det(D(l)ptD( m/ )]-iDt(i)) (3 ) 
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to represent the fermion and Pauli-Villars pieces 
of the action with a single pseudo-fermion field. 
With this approach the cancellation happens step 
by step in the leapfrog integration. We find the 
acceptance of the algorithm is increased by 10 — 
20% while the inversion costs are reduced by 20 — 
30% when using this modified force term. 

We have also implemented the chronological in- 
verter technique of 0], leading to a performance 
improvement of a factor of ks 1.7. After this im- 
provement calculating a single trajectory takes 
approximately 1.6 x 10 4 , 9 x 10 3 and 6 x 10 3 dirac 
matrix applications for m/ = 0.02, 0.03 and 0.04 
respectively. 



ffi:: 3i m m m i a 



2 4 6 



10 12 14 16 



4. RESULTS 

While the lattices collected represent part of 
a larger RBC collaboration project to calculate 
many hadronic quantities of phenomenological in- 
terest, here we will concentrate on a few mesonic 
observables to determine the basic properties of 
our simulations such as scale and quark mass. To 
calculate these quantities we have used every 50th 
trajectory, leaving out the first w 600 trajectories 
to allow the evolutions to thermalise. All quoted 
errors will be from a jackknife estimate of the sta- 
tistical error. 

To quantify the chiral symmetry breaking from 
finite L S1 we have measured the residual mass, 
defined from the breaking term in the 
Ward-Takahashi identity To extract this we 
look at 



R(t) 



(4) 



which for time greater than some t m i n should be 
time independent and equal to am res Figure 
1 shows this for the m/ = 0.02 evolution at the 
dynamical point. As can be seen, a plateau is 
evident for t>2, with an error weighted average 
between timeslice 6 and the end giving a value of 
0.00137(2). This number is relatively insensitive 
to the quark mass with linear extrapolation to 
mj = giving m res = 0.00136(5). 

While this value for the residual mass is rela- 
tively small compared to the input quark mass, 
to properly interpret the value we must know the 



Figure 1. residual mass extraction for m/ = 0.02 
evolution 
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Figure 2. M p - dynamical extrapolation 

lattice spacing and mass renormalisation. For the 
purpose of this preliminary analysis we will as- 
sign a single lattice spacing for all three evolutions 
based on a linear extrapolation in the dynamical 
quark mass of the p meson mass. This is shown 
in Fig |21 and leads to an inverse lattice spacing 
of 1.806(60)GeV. 

Figure|31shows M% for degenerate quark masses 
versus dynamical mass. Fitting to the naive, first 
order chiral perturbation theory, expectation that 



M* = B„ (mi+m 2 ) , 



(5) 



3 




5110(1 



Figure 3. M% - dynamical extrapolation Figure 4. Topological charge history for the m/ = 

0.02 evolution. 



where B n is a constant and mi and mi are 
the quark masses, gives a x 2 per degree of free- 
dom of 0.2 and results that are consistent with 
the pion mass vanishing at m/ = — m ros with 
Ml(m s = -m res ) = 1.6(11) x 10 4 MeV 2 . The 
experimental value of M.\ is shown on the figure 
as a dotted line. Together with Eq [5] this sug- 
gests our lightest quark mass is a little above half 
the strange quark mass. Going further, follow- 
ing the same approach as [7], we extract prelimi- 
nary values of the renormalised light and strange 
quark masses of 3.94(31)MeV and 103(8)MeV in 
the MS-scheme at 2 GeV. An alternative way to 
get a rough idea of the size our input quark mass 
in physical terms, which is independent of the way 
we are setting our scale, is to calculate the ratio 
Mn/Mp. This is 0.541(8), 0.598(8) and 0.637(8) 
for irif = 0.02, 0.03 and 0.04 respectively. 

In quenched simulations using the DBW2 ac- 
tion it was noticed that tunneling between differ- 
ent topological sectors is suppressed with respect 
to more standard gauge actions UJ. As such, it 
is important to study how the topological charge 
varies with trajectory number in our dynamical 
simulations. Fig 01 shows this for the mj — 0.02 
evolution, determined using a classically 0(a 4 ) 
improved definition of the topological charge cal- 
culated on each lattice after applying 20 steps of 



APE smearing with a coefficient of 0.45. While 
it is encouraging that the value of the topological 
charge is changing, it is clear that, with a sepa- 
ration of 50 trajectories between lattices, strong 
correlations are present. 

5. CONCLUSIONS 

A preliminary study of Nf = 2 dynamical 
DWF QCD, using the DBW2 gauge action at an 
inverse lattice spacing of » 1.8GeV, shows that a 
regime exists for which the explicit chiral symme- 
try breaking is small for a computationally prac- 
tical extent of the fifth dimension. 
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